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In these lectures we provide an introduction to the phase structure of QCD. We 
begin with a brief discussion of QCD, the symmetries of QCD, and what we mean 
by a "phase of QCD" . In the main part of the lectures we discuss the phase diagram 
of QCD as a function of the temperature and the baryon density. We focus, in 
particular, on the high temperature plasma phase, the low temperature and low 
density nuclear phase, and the high density color superconducting phases. 
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1. Introduction 

D ! 

(— I , In these lectures we wish to provide an introduction to the phase structure of 

QCD. The phase of QCD that we Hve in is characterized by the permanent 
confinement of quarks, and the existence of a large mass gap. There are 
rS ' several reasons for trying to understand whether other phases of QCD exist, 

' and what conditions are required in order to observe these phases: 

1) Other phases of QCD exist in the universe: The universe started out 
in a hot and dense phase. It expanded and cooled and about 10^^ sec after 
the big bang it passed through a transition from a quark gluon plasma to 
a hadronic phase. Even today, extreme conditions exist in the universe. In 
supernova explosions matter is heated to several tens of MeV, sufficient to 
evaporate ordinary nuclei. The compact remnants have central densities 
several times larger than the saturation density of nuclear matter. 

2) Exploring the entire phase diagram help us to understand the phase 
that we live in: The structure of hadrons and their interactions are deter- 
mined by the symmetries of the QCD vacuum. Studying the phase diagram 
of QCD allows us to understand the possible ways in which the symmetries 
of QCD can be realized. 

3) QCD simplifies in extreme environments: At scales relevant to 
hadrons QCD is strongly coupled and we have to rely on numerical simu- 
lations in order to test predictions of QCD. In the case of large tempera- 
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ture or large baryon density there is a large external scale in the problem. 
Asymptotic freedom implies that the bulk of the system is governed by 
weak coupling. As a result, we can study QCD matter in a regime where 
quarks and gluons are indeed the correct degrees of freedom. 

In these lectures wc will give a general introduction into the physics 
of the QCD phase diagram. There are several excellent text books and 
reviews articles that provide a much more detailed discussion of QCD and 
hadronic matter at finite temperature and density ^'^'3. We also recommend 
more specialized texts on field theory at finite temperature ^'^'^ and density 
'''^j as well as reviews on the phase structure of dense matter ^'lO'i^ and 
on color superconductivity 12,13,14,15 jj^ ^j^jg -^rite-up we will not try to 
give a summary of the experimental program at RHIC and the SPS. A 
useful reference is the series of white papers that was recently published 
by the RHIC collaborations We will also not review implications of the 
phase structure of QCD for the structure of compact stars or observational 
constraints on the behavior of dense matter io,i7,i8 

2. QCD and symmetries 
2.1. Introduction 

We begin with a brief review of QCD and the symmetries of QCD. The ele- 
mentary degrees of freedom are quark fields ip'^ ^ and gluons A'^. Here, a is 
color index that transforms in the fundamental representation for fcrmions 
and in the adjoint representation for gluons. Also, / labels the quark flavors 
u, d, s, c, 6, t. In practice, we will focus on the three light flavors up, down 
and strange. The QCD lagrangian is 

£ = ^Vi/(*^-m/)V/-^G^.G;,, (1) 
/ 

where the field strength tensor is defined by 

Gl, = d,Al-d,Al + gr'^AlAl, (2) 
and the covariant derivative acting on quark fields is 

iipi, = Y{id^+9Al^^^l^. (3) 

QCD has a number of interesting properties. Most remarkably, even though 
QCD accounts for the rich phenomenology of hadronic and nuclear physics, 
it is an essentially parameter free theory. As a first approximation, the 
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masses of the light quarks u, d, s arc too small to bo important, while the 
masses of the heavy quarks c, h, t are too heavy. If we set the masses of the 
light quarks to zero and take the masses of the heavy quarks to be infinite 
then the only parameter in the QCD lagrangian is the coupling constant, 
g. Once quantum corrections arc taken into account g becomes a function 
of the scale at which it is measured. Gross, Wilczek and Politzer showed 
that 19-20 

If the scale is large then the coupling is small, but in the infrared the 
coupling becomes large. This is the famous phenomenon of asymptotic free- 
dom. Since the coupling depends on the scale the dimensionless parameter 
g is traded for a dimensionful scale parameter Aqcd- In essence, Aqcd is 
the scale at which the theory becomes non-pcrturbative. 

Since Aqcd is the only dimensionful quantity in QCD (rriq = 0) it is 
not really a parameter of QCD, but reflects our choice of units. In standard 
units, Aqcd — 200 MeV ~ lfm~ . Note that hadrons indeed have sizes 

^QCD- 

Another important feature of the QCD lagrangian are its symmetries. 
First of all, the lagrangian is invariant under local gauge transformations 
U{x) e SU{3)c 

i){x) ^ U{x)i}{x), A^{x) ^ U{x)A^U'<{x) + iU{x)d^U\x), (5) 

where = A^(A°/2). While the gauge symmetry is intimately connected 
with the dynamics of QCD we observe that the interactions are completely 
independent of flavor. If the masses of the quarks are equal, m„ — nid = mg , 
then the theory is invariant under arbitrary flavor rotations of the quark 
fields 

V'/ ^ Vfgijg, (6) 

where V € SU{3). This is the well known flavor (isospin) symmetry of the 
strong interactions. If the quark masses are not just equal, but equal to 
zero, then the flavor symmetry is enlarged. This can be seen by defining 
left and right-handed fields 

i;L.,R^l{l±l5)i^. (7) 

In terms of L/R fields the fermionic lagrangian is 

£ = ^L{iiP)^L + ^R{iIp)il'R + tpLMiJn + ^PrM^jl, (8) 
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where M = diag(TOu, m,;, TOs). We observe that if quarks are massless, 
rriu = TTid = nrig = 0, then there is no couphng between left and right handed 
fields. As a consequence, the lagrangian is invariant under independent 
flavor transformations of the left and right handed fields. 

V'L,/ ^ LfgilL^g, i^Rj Rfgi'R,g, (9) 

where {L,R) e SU{3)l x SU{3)r. In the real world, of course, the masses 
of the up, down and strange quarks are not zero. Nevertheless, since 
'mu-i'md -^ms < Aqcd QCD has an approximate chiral symmetry. 

Finally, we observe that the QCD lagrangian has two U{1) symmetries, 

U{1)b: i^L^e'^L, i^R^e'^ijR, (10) 

U(1)a: ^L^e^^^L, V'fl^e-^«Vii- (H) 

The U{1)b symmetry is exact even if the quarks are not massless. The 
axial U{1)a symmetry is exact at the classical level but it is broken in 
the quantum theory. This phenomenon is referred to as an anomaly. The 
divergence of the U{1)a current is given by 

= ^Gl^Gl^, (12) 
where = e,ii/a/3G^^/2 is the dual fleld strength tensor. 

2.2. Phases of QCD 

The phases of QCD are related to the different ways in which the symme- 
tries of QCD can be realized in nature. We first consider the local gauge 
symmetry. There are three possible realizations of a local symmetry: 

1) Coulomb Phase: In a Coulomb phase the gauge symmetry is un- 
broken, the gauge bosons are massless and mediate long range forces. In 
particular, the potential between two heavy charges is a Coulomb potential, 
V{r) ~ 1/r. 

2) Higgs Phase: In a Higgs phase the gauge symmetry is spontaneously 
broken and the gauge bosons acquire a mass. As a consequence, the poten- 
tial between two heavy charges is a Yukawa potential, V{r) ~ exp(— rnr)/r. 
We should note that local gauge symmetry is related to the fact that we 
are using redundant variables (the four-component vector potential de- 
scribes the two polarization states of a massless vector boson), and that 
therefore a local symmetry cannot really be broken (Elitzur's theorem ^^). 
We will discuss the exact meaning of "spontaneous gauge symmetry break- 
ing" in Sect. 4.6 below. 
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3) Confinement: In a confined phase all the physical excitations are 
singlets under the gauge group. Confinement can be strictly defined only 
in theories that do not have light fields in the fundamental representation. 
In that case, confinement implies that the potential between two heavy 
charges rises linearly, V{r) ^ kr. This is called a string potential. If there 
arc light fields in the fundamental representation, as in QCD with light 
quarks, then the string can break and the potential levels off. 

It is interesting to note that all three realizations of gauge symmetry 
play a role in the standard model. The U{1) of electromagnetism is in a 
Coulomb phase, the SU{2) is realized in a Higgs phase, and the SU{3) of 
color is confined. Also, as we shall see in these lectures, there are phases of 
QCD in which the color symmetry is not confined but realized in a Higgs 
or Coulomb phase. 

Different phases of matter, like liquid vs solid, superfluid vs normal, are 
related to the realization of global symmetries. In QCD at zero baryon 
density spacetime symmetries as well as U{1) symmetries cannot be broken 
22,23 This means that phases of QCD matter are governed by the real- 
ization of the chiral SU{3)l x 5f7(3)i{ symmetry. If the baryon density is 
not zero both space-time and U{1) symmetries can break and the phase 
diagram is much richer. 

2.3. The QCD vacuum 

In the QCD ground state at zero temperature and density chiral symmetry 
is spontaneously broken by a quark- anti-quark condensate {'tptp)- We can 
view the chiral condensate as a matrix in flavor space. In the QCD vacuum 

{^i^R) = {i^yi) - -5^^(230 MeV)^ (13) 

which implies that chiral symmetry is spontaneously broken according to 
SU{3)l X SU{i)R SU{3)v The SU{3)v flavor symmetry is broken 
explicitly by the difference between the masses of the up, down and strange 
quark. Since rus ^ mu,md the SU{2) isospin symmetry is a much better 
symmetry than SU (3) flavor symmetry. 

Chiral symmetry breaking has important consequences for the dynamics 
of QCD at low energy. Goldstone's theorem implies that the breaking 
of SU{3)l X SU{3)r SU{3)v is associated with the appearance of an 
octet of (approximately) massless pseudoscalar Goldstone bosons. Chiral 
symmetry places important restrictions on the interaction of the Goldstone 
bosons. These constraints are most easily obtained from the low energy 
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effective chiral lagrangian. The transformations properties of the chiral 
field E follow from the structure of the chiral order parameter, 

S^LSijt, St^iJE^it, (14) 

for {L,R) e SU{3)l x SU{3)r. In the vacuum we can take (S) = 1. 

Goldstone modes are fluctuations of the order parameter in the coset space 
SU{3)l X SU{3)bJ SU{3)v- They are parameterized by unitary matrices 
S = exp(2A°0"//^) where A" are the Gell-Mann matrices and = 93 
MeV is the pion decay constant. At low energy the effective lagrangian 
for S can be organized as an expansion in the number of derivatives. At 
leading order in (cJ/Ztt) there is only one structure which is consistent with 
chiral symmetry, Lorentz invariance and C,P,T. This is the lagrangian of 
the non-linear sigma model 

£ = :^Tr [a^Sa^St] + . . . . (15) 

In order to show that the parameter is related to the pion decay am- 
plitude we have to gauge the non-linear sigma model. This is achieved by 
introducing the gauge covariant derivative V^S = d^'E + i^^W^S where 
W/j, is the charged weak gauge boson and is the weak coupling con- 
stant. The gauged non-linear sigma model gives a pion-V7 boson interac- 
tion C = gwfTrWj^d^TT^ which agrees with the standard definition of in 
terms of the pion- weak axial current matrix element. 

Expanding S in powers of the pion, kaon and eta fields cp'^ we can derive 
low energy predictions for Goldstone boson scattering. In the pion sector 
we have 

^ = \id,rr + ^ [ird.rr - irno.^'r] +0(71)' (16) 

which shows that the low energy TTTT-scattering amplitude is completely 
determined by f,^ . Higher order corrections originate from loops and higher 
order terms in the effective lagrangian. 

In QCD chiral symmetry is explicitly broken by the quark mass term 
ij)Mip, where M = diag(TO„, m^, m^) is the quark mass matrix. In order 
to determine how the quark masses appear in the effective lagrangian it 
is useful to promote the mass matrix to a field which transforms as M — > 
LMR'^ under chiral transformations. This means that the mass term is 
chirally invariant and explicit breaking only appears when M is replaced 
by its vacuum value. There is a unique term in the chiral lagrangian which 
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is SU{3)l X SU{3)r invariant and linear in M. To order 0{d'^,M) the 
effective lagrangian is 



The mass term acts a potential for the chiral field. We observe that if the 
quark masses arc real and positive then the minimum of the potential is 
at (E) = 1, as expected. If some of the quark masses are negative unusual 
phases of QCD can appear, see 

The vacuum energy is Eyac = — 2BTr[M]. Using (V'V') = dEyac/{dm) 
wc find (ijjijj) = —2B. Fluctuations around the vacuum value S = 1 deter- 
mine the Goldstone boson masses. The pion mass satisfies the Gell-Mann- 
Oaks-Renner relation 



and analogous relations exist for the kaon and eta masses. 
2.4. QCD vacuum for different Nc and Nf 

QCD is a strongly interacting gauge theory with almost massless quarks. 

It seems natural that in such a theory bound states of quarks and anti- 
quarks are formed, that bound states in the scalar channel condense, and 
that chiral symmetry is broken. But even if chiral symmetry breaking is not 
surprising, it is not a priori clear whether the observed pattern of chiral sym- 
metry breaking and confinement is required on general grounds, or whether 
it is a particular dynamical feature of QCD. Some obvious questions are: 
Are all asymptotically free gauge theories confining? Does confinement im- 
ply chiral symmetry breaking (or vice versa)? Is the symmetry breaking 
pattern SU{3)l x SU{S)r SU{3)v unique? 

An interesting context in which these questions can be studied is the 
phase diagram of QCD and supersymmetric generalizations of QCD as a 
function of Nc and Nj, see Fig. 1. For our purposes supersymmetric QCD is 
simply a QCD-like theory with extra fermions in the adjoint representation 
and extra colored scalar fields. Including supersymmetric theories is useful 
because supcrsymmetry provides additional constraints that determine the 
symmetries of the ground state. The following interesting results have been 
obtained: 

1) In supersymmetric QCD there is a window A^c + 1 < A/ < 3A^c in 
which the theory is asymptotically free but not confining There are 
several reasons to believe that such a window exists in QCD, too. One is 




(17) 



(18) 
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Figure 1. Ground state of QCD and SUSY QCD as a function of Nc and Nf. The 
symmetry breaking pattern in SUSY QCD was clarified in a series of papers by Seiberg 
and collaborators. The phase structure of QCD is an educated guess, discussed in more 
detail in the review (Schafer and Shuryak, 1998). 



the fact that as a function of the number of flavors the second coefficient of 
the beta function changes sign before the first one does In this regime 
the coupUng constant flows to a finite value at large distance and the theory 
is scale invariant. 

2) Supersymmetric QCD also provides examples for theories that have 
confinement but no chiral symmetry breaking. This happens for Nj = 
Nc + 1. This theory contains both massless mesons and massless baryons. 
An important constraint is provided by the 't Hooft anomaly matching 
conditions 27,28^ QCD these relations show that confinement without 
chiral symmetry breaking is a possibility for Nf — 2, but ruled out for 
Nf>2. 

3) The 't Hooft consistency conditions also provide constraints on the 
symmetry breaking pattern. In QCD these conditions are not sufficiently 
strong to fix the ground state completely, but one can show that SU(3)l x 
SU{3)r SU{3)v is favored in the limit N^ ^ oo ^9. 
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4) One can show that in QCD chiral symmetry breaking imphcs a non- 
zero quark condensate In particular, one can rule out the possibility 
that (^^V) = 0, but {{'4jtpf) + 0. 

3. QCD at finite Temperature 
3.1. General Arguments 

In this section we shall discuss the phase structure of QCD at non-zero 

temperature. We begin by reviewing some general arguments in favor of 
the existence of a critical temperature Tc above which quarks and gluons 
are deconfined and chiral symmetry is restored. 

Asymptotic freedom clearly suggests that the high temperature phase is 
a weakly interacting plasma Consider a non-interacting gas of quarks 

and gluons at high temperature. The typical momenta are on the order of 
the temperature, p 3T, and the density is n ~ T^. Now imagine that we 
turn on the coupling. Does this lead to a qualitative change in the system? 
Quarks and gluon can scatter but since the typical momenta are large a 
significant change in the momentum of the scattered particles requires a 
large momentum transfer. Asymptotic freedom implies that the effective 
coupling at this scale is small, and that large angle scattering events are 
rare. If the change of momentum is small then the scattering involves 
large distances and the interaction is modified by the dense medium. We 
will see below that the quark-gluon medium screens the interaction and 
that the effective interaction is again weak. There is a small subtlety here, 
as static magnetic interactions are not screened. This implies that high 
temperature QCD has a genuinely non-perturbative sector, but this sector 
is not important as far as bulk properties of the high temperature phase are 
concerned. We conclude that the assumption of a weakly interacting quark- 
gluon system at high temperature leads to a self consistent picture. Since 
this system will exhibit medium effects, such as damping and screening, 
collective modes, etc. that are typical of plasmas it was termed the quark 
gluon plasma (QGP) 

It is instructive to consider a simple model of the equation of state. The 
pressure and energy density of non-interacting massless particles is 



where g is the number of degrees of freedoms. In a quark gluon plasma we 
have gq = ANfNc quarks and gg = 2{N^ — 1) gluon degrees of freedom. For 




(19) 
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Nf = 2 we get geff = gg + = 37 and 

P = (20) 
90 ^ ' 

At low temperature the relevant degrees are Goldstone bosons. Near Tc 

we can assume that Goldstone bosons are approximately massless and the 

number of degrees of freedom is g = {Nj — 1). For Nf = 2 we get g = 3 

and 

P=— (21) 

This result seems to show that the pressure in the low temperature phase 
is always smaller than the pressure in the high temperature phase. This 
cannot be right as it would imply that the phase with chiral symmetry 
breaking is never favored. The problem is that there are non-perturbative 
effects in the low temperature phase. These effects give rise to a negative 
vacuum energy and a positive vacuum pressure. Lorentz invariance implies 
that the vacuum energy momentum tensor is of the form T^j, = Bgf^^ and 

^vac — Pvac — ~\~B. (^^) 

In QCD the vacuum energy satisfies the trace anomaly relation 

evac = -^{-G^} ^ -0.5 GeV/fm^. (23) 

The numerical value comes from QCD sum rule determinations of the gluon 
condensate {ctG"^) and has a considerable uncertainty. We can now obtain 
an estimate of the transition temperature by requiring that the pressure in 
the low temperature phase equals the pressure in the quark gluon phase. 
We find 

= ( 17^ j ^ (^^^ 

We can also determine the critical energy density. The energy densities just 
below and just above the critical temperature are given by 

e(T-) = —T^ ^ 130 MeV/fm^ (25) 

e(r+) = —-T^ + 3-:^ 2000 MeV/fm^ (26) 
ou 

We observe that the energy density in the QGP exceeds 1 GeV/fm'^. This 
should be compared to the energy density in cold nuclear matter which is 
about 150 MeV/fm^ 
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An independent estimate of the transition temperature c;an be obtained 
using the chiral effective theory. We saw that the chiral condensate is 
related to the mass dependence of the vacuum energy. At tree level and to 
leading order in the quark masses the condensate is given by the coefficient 
B in the chiral lagrangian. We can also calculate corrections to this result 
due to thermal Goldstone boson fluctuations. At leading order it is sufficient 
to consider the free energy of a non-interacting pion gas 

where E^^ = y/p^ + m^. The quark condensate is {'(pip) = {Nf)~^dF/dm. 
Equation (27) depends on the quark mass only through the pion mass. 
Using the Gell-Mann-Oakes-Renncr relation (18) we find 

{^h = mo |l - {^^) + . . j . (28) 

This result indicates that the chiral condensate vanishes at a critical tem- 
perature 



- - 200 MeV {Nf = 3), (29) 

which is roughly consistent with the estimate obtained in equ. (24). 

3.2. Chiral Symmetry Restoration 

In the vicinity of the phase transition QCD is genuinely non-perturbative 

and it is hard to improve on the rough estimates provided in the previous 
section. One possibility for determining the transition temperature and 
elucidating the nature of the phase transition is the use of large scale nu- 
merical simulations. We will discuss this approach in Sect. 3.4. Before we 
do so we would like to review certain general statements about the phase 
transition that follow from the symmetries of the low and high temperature 
phase. 

We begin with the chiral phase transition. We shall assume that the 
chiral transition is a second order phase transition, i.e. the order param- 
eter goes to zero continuously. We will explore the consequences of this 
assumption and check whether it leads to a consistent picture. In the vicin- 
ity of a second order phase transition the order parameter fluctuates on 
all length scales and the correlation length diverges. This means that the 
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Table 1. Correspondence between the chiral phase transition in QCD 
and the ferromagnetic transition in a four-component magnet. 

SU{2)l X SU{2)ii QCD 0(4) magnet 

{■4>ij}) X condensate M magnetization 

niq quark mass H3 magnetic field 

TT pions </> spin waves 

details of the interaction are not important, and only the symmetry of the 
order parameter matters. 

In QCD with two flavors the order parameter is a 2 x 2 matrix U^^ = 
(ip-^ipfi) . We can define a four component vector cf)'^ by U^^ = (/)°(r°)-'^ with 
r" = (f, 1). Chiral transformations {L,R) S SU{3)l x SU{3)ii correspond 
to rotations (j)" — > RF'^'cf)^ with R G 50(4). In the low temperature phase 
chiral symmetry is broken and (0") = cr6'^'^. Near the order parameter is 
small and we we can expand the free energy in powers of the order parameter 
and its derivatives. Chiral symmetry implies that only the length of can 
enter. To order and to leading order in gradients of the fields we have 

F = y'ci'^{^(vrr + Y(W + ^(W + •••}, (30) 

where /i^,A are parameters that depend on the temperature. Equ. (30) 
is the Landau-Ginzburg effective action. Note that fluctuations of the or- 
der parameter are dominated by static fields (f)°'{x,t) = <j)°'{x). This will 
be explained in more detail in Sect. 3.5. The main observation is that 
fiuctuations with energy much smaller than ttT are described by a three 
dimensional theory. 

Stability requires that A > 0. Below the critical temperature /x^ < 
and chiral symmetry is broken. At Tc the parameter fj,^ changes sign and 
we can write /z^ = n^t where t = (T — T^jT^ is the reduced tempera- 
ture. As a first approximation we can ignore fiuctuations and study the 
Landau-Ginzburg action in the mean field approximation. In that case the 
chiral order parameter goes to zero as {<^) ^ t^/'^. This result is modi- 
fied if fluctuations are included. This can be done using rcnormalization 
group methods or numerical simulations. These methods also demonstrate 
that near higher order operators not included in equ. (30) are indeed 
irrelevant. The results are 

C-i-° a = -0.19, 

(V^V) P = 0.38, (31) 

m^r^t" V = 0.73, 
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where C is the specific heat and the Goldstonc boson mass is defined 
as the inverse correlation length of spatial fluctuations of the pion fleld. 

The coefficients a, /3, v are called the critical indices of the phase transi- 
tion. These coefficients are universal, i.e. independent of the details of the 
microscopic details. For example, the critical indices of the chiral phase 
transition in QCD agree with the critical indices of a four-component mag- 
net in d = 3 space dimensions, see Table 1. 

In QCD with Nf = 3 flavors the order parameter is a 3 x 3 matrix 
U^^ . The main new ingredient in the Landau-Ginzburg theory is a chirally 
invariant cubic term det(t/) + h.c. It is easy to check that the cubic term 
will lead to an effective potential that has two degenerate minima at the 
transition temperature. This implies that the transition is first order and 
the order parameter goes to zero discontinuously. Even if the coefficient 
of the cubic term is zero initially, fluctuations tend to generate a cubic 
interaction as T ^ Tc- A transition of this type is called fluctuation induced 
first order. 



3.3. Deconfinement 

It is not immediately obvious how to identify the symmetry associated with 
the deconfinement transition. In order to make contact with the lattice for- 
mulation discussed in Sect. 3.4 wc will study this question in euclidean 
space, i. e. with the time variable analytically continued to imaginary time, 
t it. We first derive an expression for the potential between two heavy 
quarks. Consider the Dirac equation for a massive quark in a gluon back- 
ground field 

(do - igAo - a{iV + gl) + 7oM) V = 0. (32) 

In the limit mg — > oo we can ignore the spatial components of the covariant 
derivative. The quark propagator is 

S{x,x') ^ exp (^ig J Aod?j e-™(*-*')j(f - x'). (33) 

The heavy quark potential is related to the amplitude for creating a heavy 
quark-anti-quark pair at time t = 0, separating the pair by a distance R, 
and finally annihilating the two quarks at time t = T. The amplitude is 
related to the Wilson loop 

W{R, T) = exp [ig <f A^^dzA , (34) 
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where the integration contour is a i? x T rectangle and wc have dropped the 
Dirac projection operators (l+7o)/2. If T i? the ampUtude is dominated 
by the ground state and we expect W{R,T) ~ exp{—V{R)T) where V{R) 
is the heavy quark potential. Confinement implies that V{R) ~ kR and 
the Wilson loop satisfies an area law 

W{R,T)r^exp{-kA),. (35) 

where A = RT. In order to construct a local order parameter we consider 
the Polyakov line 

P{x) = -^^T^im] = ^^'Tr 

As we will explain in more detail in the next section gauge fields are pe- 
riodic in imaginary time with period /3 = 1/T. The Polyakov line can be 
interpreted as the free energy of a single quark. (P) ~ exp(— mg/?). In the 
confined phase we expect that the energy of an isolated quark is infinite, 
while in the deconfined phase it is finite. This implies that 

(P) = confined, (P) ^ deconfined. (37) 

The global symmetry of the order parameter is P — > zP with z = 
exp{2'jTki/Nc) € Zn^. Since Zn^ is the center of the gauge group SU{Nc) 
this is sometimes called the center symmetry. Color singlet correlation 
functions always involve combinations of Polyakov lines that are invariant 
under center transformations. A heavy baryon correlation function, for 
example, is of the form tr[{L{x))'^^] and is invariant because z^" = 1. A 
non-zero expectation value for the Polyakov line, on the other hand, breaks 
center symmetry. 

Note that the symmetry is broken in the high temperature phase and 

restored in the low temperature phase. This might seem somewhat unusual, 
but there are examples of spin systems that have an equivalent "dual" 
description in terms of a gauge theory In the dual theory the high and 
low tempcratiire phases are interchanged. 

The Zn^ Landau-Ginzburg action is given by 

F = yrf-\- |ilVPp + /i2|P|2+5Re(p3) + A|P|* + ...|. (38) 

The cubic term is allowed only if TVg = 3. As in the case of the chiral phase 
transition we expect the cubic term to drive the transition first order. The 
two color theory is in the equivalence class of the Z2 Ising model, which is 
known to have a second order transition. The three color theory is in the 



exp ig 



Andt 



(36) 
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equivalence class of a three state Potts model, which does indeed show a 
first order transition. 

The phase structure as a function of the light quark mass m = m„ = rn^ 
and the strange quark mass rris is summarized in Fig. 2. The lower left hand 
corner of the diagram is m = = and corresponds to three massless 
quarks. In that case we expect a first order chiral phase transition. Along 
the diagonal we have m — and the SU{3)v flavor symmetry is exact. As 
the quark masses increase the strength of the first order transition becomes 
weaker and the transition eventually ends at a second order critical point. 
If the light quarks are kept massless as the strange quark mass is increased 
the endpoint of the first order transition is a tricritical point at which the 
the endpoint of the first order transition in the three flavor theory meets 
the second order transition of the two flavor theory. This transition turns 
into a smooth crossover as soon as the light quarks are given a small mass. 

The upper right hand corner of the plot is m = rus — > oo and corre- 
sponds to the pure glue theory which has a first order transition. Dynamical 
quarks break the Z^^ symmetry and the strength of the first order transi- 
tion decreases as the quark masses are lowered. The first order transition 
in the pure glue theory ends on a line of second order transitions. 

We do not know with certainty where the physical point is located on 
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Figure 3. Chiral and deconfinement transitions for Nf = 2 dynamical quark flavors, 
from (Karsch 2002). The two figure show the chiral condensate and the Polyakov line 
as a function of the bare coupling /3 = (In this figure, /3 is not 1/T). Asymptotic 

freedom implies that on a fixed lattice Nt x A^^ increasing /3 corresponds to increasing 
the temperature. The susceptibilities Xrm XL a-re related to order parameter fluctuations, 
e.g. XL=nU{L^)-{m. 



this phase diagram. Lattice calculations currently favor the possibility that 
the phase transition is in the crossover region, closer to the first order chiral 
transition than to the first order deconfinement transition. 

We should emphasize that Fig. 2 focuses on regions of the phase diagram 
in which there is a sharp phase transition and the order parameter is a non- 
analytic function at Tc- The figure should not be taken to imply that the 
chiral and deconfinement transitions are completely separate phenomena, 
or that the transition cannot be observed in the crossover region. Fig. 3 
shows the order parameter as well as the order parameter susceptibilities 
for the chiral and deconfinement transitions. The results were obtained 
from lattice calculations with semi-realistic values of the quark masses. 
We observe that even though both transitions are crossovers clear peaks 
in the susceptibilities are still visible. We also note that the chiral and 
deconfinement transitions occur at the same temperature. This can be 
understood in models that include a non-zero coupling between the effective 
actions for the chiral and deconfinement order parameters 39,40,4i^ 



3.4. Lattice QCD 

Symmetry arguments cannot determine the critical temperature, the crit- 
ical energy density, and many other properties of matter near the phase 
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transition. In order to compute these properties we have to rely on numer- 
ical simulations of the QCD partition function 

Z = Tj:[e-'^"], P = l/T F = TV-llog{Z), (39) 

where H is the QCD Hamiltonian, f3 is the inverse temperature and F is the 
free energy. We can write the partition function as a quantum mechanical 
evolution operator Z = Tr[e~'^~*''^^] with imaginary time r = — The 
evolution operator can be written as a path integral 




where £b is the imaginary time (euclidean) lagrangian and we have to 
impose (anti)periodic boundary conditions on the quark and gluon fields 

A^{x,p)=A^{x,0), V(^,/J) = -V(^,0). (41) 

The path integral is an infinite dimensional integral. In order to perform 

numerical simulations wc have to discrctizc space and time and introduce 
a iVx X lattice with lattice spacing a. In order to maintain exact gauge 
invariance on a discrete lattice the gauge fields are discretized in terms of 
the link variables 

?7^(n) = exp{igaA^{n)), (42) 

where n = [rir, nx,ny,nz) labels lattice sites and /U = 1, . . . , 4. In terms of 
the link variables it is easy to define a gauge covariant discrete derivative 

i[t/_4n)(/.(n + M)-</.(n)], (43) 

where n + n is the lattice site reached by starting at n and doing one hop 
in the n direction. (j){n) is a scalar field in the fundamental representation. 
The action of the pure gauge theory is given by 

5 = 4 I^Re Tr [1 - [/^(n)t/,(n + n)U-^{n + + J^)U-^{n + p)] , (44) 

which involves a loop of link variables called a plaquette. It is easy to 
check that equ. (44) reduces to the continuum action as a — * 0. Fermion 
fields tp{n) are discretized on lattice sites. There are some subtleties in 
discretizing chiral fermions which we do not discuss here The action is 
bilinear in the fermion fields. This means that the integral over the fermion 
fields can be done exactly and we are left with the determinant of a matrix 
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that depends only on the Unk variables. The lattice representation of the 
partition function is 

Z = fY[dU^{n)det{M{U))e-^, (45) 

where M{U) is the fermion matrix. The partition function depends on 
the number of lattice sites Nt-jN^-, the bare coupling constant <?, and the 
dimensionless quark masses ma. 

Note that the partition function has no explicit dependence on the lat- 
tice spacing. Asymptotic freedom implies that the bare coupling should go 
to zero as a — > and the continuum limit corresponds to 5 — > 0. Asymp- 
totically 

aKiat = exp(-87rV(65')), (46) 

where b is the first coefficient of the beta function, see equ. (4), and Aiat 
is the QCD scale parameter on the lattice. Ajat can be related to the 
continuum scale parameter by a perturbativc calculation . In practice 
the lattice spacing is not small to enough to use the perturbativc result 



February 2, 2008 3:55 Proceedings Trim Size: 9in x 6in 



hugsproc 



19 




P 



Figure 5. One loop contribution to the photon polarization tensor (left panel) and 
plasmon dispersion relation in a hot QED plasma. 

equ. (46) and a is determined by from a physical quantity like the string 
tension or the rho meson mass. Once the lattice spacing is known the 
temperature is determined by T = l/{Nra). 

Lattice results for the order parameter and the equation of state are 
shown in Figs. (3,4). Current results for the transition temperature are 



Tc(A^/ = 2) = (173±8) MeV, Tc(7V/ = 0) = (271 ± 2) MeV, (47) 



where the errors are purely statistical. The equation of state shows a rapid 
rise in the pressure near Tc, but the pressure remains significantly below 
the free gas limit even at T ~ (2 — 3)T'c. 

3.5. Perturbative QCD 

At temperatures significantly above Tc quarks and gluons are weakly cou- 
pled and perturbative methods are expected to be useful. The starting 
point of the perturbative expansion at non-zero temperature is the path in- 
tegral representation of the QCD partition function given in equ. (40) . The 
only difference as compared to the zero temperature result is the fact that 
the fields satisfy periodic boundary conditions in imaginary time. Consider 
the Fourier representation of the gluon field 



We can also write down an analogous expansion for fermions. The bound- 
ary conditions imply that the allowed frequencies, called the Matsubara 




n 



(48) 
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frequencies, are discrete. We have 

cOn = 2imT bosons , , 

ujn = (2n + l)irT fermions ' 

The only change in the Feynman rules is that continuous energy variables 
are replaced by discrete Matsubara frequencies, po — > and that inte- 
grals over energy are replaced by discrete sums 

d^p f cPp 



J (27r)4^ 4-7 (27r)3- 



(50) 



Typical sums that appear in one-loop calculations in the Matsubara for- 
malism are 

1 _ 27r /I 1 \ 

^ x2 + A;2 ~ T U ^ e2-- - ij ' ^^^^ 

k 

^ a;2 + (2fc+l)2 " X (2 " e'^^ + l) ' ^^^^ 

We observe that performing sums over Matsubara frequencies leads to Bose- 
Einstein and Fermi-Dirac distribution functions. 

As an application of the finite temperature formalism we wish to study 
the one-loop correction to the gluon propagator in a hot QCD medium. For 
simplicity we begin with the analogous problem in a hot QED plasma. The 
photon polarization function is (see Fig. 5) 

^tr[7^A7.(^ - ^)]A(fc)A(A: - q), (53) 

where A{k) = u^ + k^. Using identities like equ. (52) we can decompose the 
integral into a T = and a finite temperature part. In the following we will 

only consider the T 7^ terms. Thermal corrections to the photon propa- 
gator become important when the photon momentum is much smaller than 
the temperature. In this case we can assume that the loop momenta are on 
the order of T. This is called the hard thermal loop (HTL) approximation 
The photon polarization function in the HTL approximation is 

n,. = 2m^/^(^^^ + <5,o^, (54) 

where = e^T^/B, K = (— i, k) and is an integral over the direction 
of k. In the case of the gluon propagator in QCD there are extra graphs 
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Figure 6. One loop contribution to the gluon polarization tensor in the quark gluon 
plasma. Solid lines are quark propagators, wavy lines are gluons, and dashed lines are 
ghosts. 



generated by the three and four-gluon vertices as weU as ghost contribu- 
tions but, remarkably, the structure of the HTL polarization function is 
unchanged. In QCD the parameter is given by = g^T^{l + Nf /6). 

Insertions of the polarization function into the photon (gluon) propa- 
gator form a simple geometric series. The resummed photon propagator 
is 

This result can be used to study the interaction between two charges in the 
plasma. The Coulomb interaction is determined by the Fourier transform 
of the static propagator 



(27r)'^ + IIoo r 



where mjj = 2m^ is the Debye mass and we have used noo(0, 0) = mjj. 
Equ. (56) shows that the Coulomb interaction is screened at distances r ~ 
m^)^ ~ (eT)^^. The mechanism for charge screening is easy to understand. 
A test charge polarizes the electron-positron plasma, and the polarization 
cloud screens the charge. Note that the typical distance between charges 
is r ~ ^ ru and Debye screening is a collective effect that involves 
many charge carriers. 

The magnetic interaction is not screened, 11^(0, (f ^ 0) = 0. However, 
if the photon energy is finite the polarization develops an imaginary part 

lmnii{uj,q) —mjjQiq - u), (57) 

q 

and non-static exchanges are damped. This phenomenon is called Landau 
damping. The mechanism of Landau damping is a transfer of energy from 
the electromagnetic field to electrons and positrons in the plasma. The 
absence of magnetic screening implies that the static magnetic sector of 
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the QCD plasma remains non-perturbative even if the temperature is very 
high. 

In order to study the propagation of cohective modes in the plasma in 
more detail it is useful to split the polarization tensor into transverse and 
longitudinal components 

n^(g)Pj,+n^(g)P^^, (58) 
Sij-Mj, P^o = Poi=0, (59) 

-5.. + ^ - Pj.- (60) 

We can study the propagation of photons by identifying the poles of the 
transverse and longitudinal components of the photon propagator, see 
Fig. 5. Wc observe that for large momenta \q\ ^ m the dispersion re- 
lation is not strongly affected by the medium. In this limit we also find 
that the longitudinal mode has an exponentially small residue. As ^ ^ 
the energy of both longitudinal and transverse modes approach the plasma 
frequency cOpi = -\/2/3 m. 

4. QCD at Small Density: Nuclear Matter 
4.1. Introduction 

In this section we study hadronic matter at non-zero baryon density. In 
QCD the numbers of all quark flavors are conserved individually. Once 
the weak interaction is taken into account only baryon number and electric 

charge arc conserved. Bulk matter has to be electrically neutral because the 
Coulomb energy of a charged system diverges in the infinite volume limit. In 
hadronic matter neutrality can be achieved by balancing the charge density 
in hadrons, which is usually positive, by a finite density of electrons. 

The partition function of QCD at non-zero baryon chemical potential 
is given by 

i ^ ^ 

where i labels all quantum states of the system, Ei and iVj are the energy 
and baryon number of the state i. If the temperature and chemical poten- 
tial are both zero then only the ground state contributes to the partition 
function. All other states give exponentially small contributions. QCD has 
a massgap for states with non-zero baryon number. This means that there 



n^.(g) 
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Figure 7. Naive phase diagram of hadronic matter as a function of the baryon chemical 
potential and temperature. 

is an onset chemical potential 

^lonset = niin(£'j/iVj), (62) 

i 

such that the partition function is independent of /i for \x < fionset ■ For fi > 
fJ-onset the baryon density is non-zero. If the chemical potential is just above 
the onset chemical potential we can describe QCD, to first approximation, 
as a dilute gas of non-interacting nucleons. In this approximation [lonset — 
■niN- Of course, the interaction between nucleons is essential. Without 
it, we would not have stable nuclei. As a consequence, nuclear matter is 
self-bound and the energy per baryon in the ground state is given by 

E 

-mAT ~ -15MeV. (63) 

The onset transition is a first order transition at which the baryon density 
jumps from zero to nuclear matter saturation density, po — 0.14fm^''. The 
first order transition continues into the finite temperature plane and ends 
at a critical endpoint at T = Tc ~ 10 MeV, see Fig. 7. 

Nuclear matter is a complicated many-body system and, unlike the situ- 
ation at zero density and finite temperature, there is little information from 
numerical simulations on the lattice. This is related to the so-called 'sign 
problem'. At non-zero chemical potential the path integral representation 
of the partition fmiction is 



det(i^ + i^74)e"'^ = j dA^e"l'\Aei{iIp + i^M"i4)\e~^ , (64) 
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where (j) is the complex phase of the fermion determinant. Since the deter- 
minant is complex standard Monte-Carlo techniques based on importance 
sampling fail. Recently, some progress has been made in simulating QCD 
for small ^, and T ~ Tc ^5,46,47^ ^j^g regime of small temperature remains 
inaccessible. 

However, if the density is very much larger than nuclear matter satu- 
ration density, p ^ po, we expect the problem to simplify. In this regime 
it is natural to use a system of non-interacting quarks as a starting point 
The low energy degrees of freedom are quark excitations and holes 
in the vicinity of the Fermi surface. Since the Fermi momentum is large, 
asymptotic freedom implies that the interaction between quasi-particles is 
weak. As a consequence, the naive expectation is that chiral symmetry is 
restored and quarks and gluons are dcconfined. It seems natural to assume 
that the quark liquid at high baryon density is continuously connected to 
the quark-gluon plasma at high temperature. These naive expectations are 
summarized in the phase diagram shown in Fig. 7. 

4.2. Fermi liquids 

Before we study the high density phase in more detail we would like to 
discuss systems of nucleons at low density. For simplicity we will begin with 
pure neutron matter at densities below nuclear matter saturation density. 
This problem is relevant to the behavior of matter near the surface of a 
neutron star, which is at subnuclear densities and has a large neutron-to- 
proton ratio. We will also see that pure neutron matter exhibits some very 
interesting universal features which can be studied experimentally using 
trapped atomic gases. 

If the density is low then the typical momenta are small and neither the 
structure of the neutron nor the details of the neutron-neutron interaction 
are important. This means that the system can be described by an effective 
lagrangian of pointlike nucleons interacting via a short-range interaction 
The lagrangian is 



The coupling constant Co is related to the scattering length, Co = ^ira/m. 
Note that Co > corresponds to a repulsive interaction, and Co < is an 
attractive interaction. The lagrangian equ. (65) is invariant under the U (1) 
transformation if) e^'^ip. The U{1) symmetry implies that the fermion 




(65) 
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Figure 8. Leading order Feynman diagrams for the ground state energy of a dilute gas 
of fermions interacting via a short range potential. 



number 



N= d?x^p^ 



(66) 



is conserved. We introduce a chemical potential /i conjugate to the fermion 
number N and study the partition function 



=Tr 



(67) 



Here, H is the Hamiltonian associated with C and (3 = 1/T is the inverse 
temperature. The average number of particles for a given chemical potential 
/i and temperature T is given by (A'') = T{d\ogZ)/{dy). At zero temper- 
ature the chemical potential is the energy required to add one particle to 
the system. 

We observe that the chemical potential simply shifts the energy in the 
lagrangian. This implies that we have to carefully analyze the boundary 
conditions in the path integral in order to fix the pole prescription. The 
correct Minkowski space propagator is 



where e„ = 



= Saf 



Pq 



fi, Ep = p^/(2to) and S 



i6 
0+. 



Pq 



i6 



(68) 



The quantity pp — 
y^2mfi is called the Fermi momentum. The two terms in equ. (68) have 
a simple physical interpretation. At finite density and zero temperature 
all states with momenta below the Fermi momentum are occupied, while 
all states above the Fermi momentum are empty. The possible excitation 
of the system are particles above the Fermi surface or holes below the 
Fermi surface, corresponding to the first and second term in equ. (68). The 
particle density is given by 



N 
V 



d'^p 



+0+ 



d^p 
(2^ 



QiPF -p) 



37r2' 



(69) 
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Figure 9. Particle-particle ladder diagrams for the ground state energy of a dilute gas 
of fermions. 



As a first simple application wc can compute the energy density as a func- 
tion of the fermion density. For free fermions, we find 

We can also compute the corrections to the ground state energy due to 
the interaction ^Co{'tp^tp)'^. The first term is a two- loop diagram with one 
insertion of Co, see Fig. 8. We have 

A=C.(ii)\ (71) 

We should note that equ. (71) contains two possible contractions, called the 
direct and the exchange term. If the fermions have spin s and degeneracy 

g = (2s + 1) then cqu. (71) has to be multiplied by a factor g(g — 1) /2. Wc 
also note that the sum of the first two terms in the energy density can be 
written as 

E Pf 2 



which shows that the Co term is the first term in an expansion in ppo,, 
suitable for a dilute, weakly interacting, Fermi gas. The expansion in (ppa) 
was carried out to order {ppa)'^ by Huang, Lee and Yang ^9,50 gingg then, 
the accuracy was pushed to '"'^^ 0((pira)*log(pira)). 



4.3. Unitary Limit 

The neutron- neutron scattering length is very large, a„„ = — 18 fm, and 
the (ppa) expansion is not particularly useful. Indeed, since the scattering 
length is so much larger than all other hadronic length scales it makes sense 
to consider the opposite limit and take the scattering length to infinity. This 
means that there is a bound state right at threshold and that the low energy 
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Figure 10. Total energy of an interacting fcrmion gas in units of the energy of a free 
fermion gas as a function of {kpa)~^. The open circles show the result of a numerical 
calculation of particle-particle ladder diagrams. The dashed curve shows the approxima- 
tions given in equ. (74). The star is the result of the d — > oo calculation in the unitary 
limit. 



cross section saturates the unitarity bound. If neutron matter is dilute then 
we can also assume that (pft) <C 1, where r is the range of the potential. 
In this limit the only energy scale in the problem is p|n/(2m) and the energy 
per particle is 

^-fM (73) 

where ^ is an unknown parameter. Comparison with equ. (70) shows that 
for free fermions ^ = 1. 

Neutron matter in the unitary limit is very strongly correlated and the 
determination of f is a complicated, non-perturbative problem. However, 
since ^ is insensitive to the details of the interaction the result is the same 
for any dilute Fermi gas with a two-body bound state near threshold. It 
is now possible to create such a system in the laboratory by trapping cold 
fermionic atoms. In these systems the scattering length can be controlled 
using Feshbach resonances induced by an external magnetic field. A recent 
experimental analysis yields the value ~ 0.45. 

There have been a number of theoretical attempts to determine ^. Since 
the two-body interaction is large it is natural to begin with the sum of all 



Februciry 2, 2008 3:55 Proceedings Trim Size: 9in x 6in 



hugs'proc 



28 

two-body diagrams, see Fig. 9. This sum gives 

^ _ pj. [3 2(fcfa)/(37r) 
N 2M\5 l-^{n-2\og{2)){ppa) 

from which we deduce ^ ~ 0.32. This is reasonably close to the experimen- 
tal result, but since the system is strongly correlated there is no obvious 
reason to restrict ourselves to two-body ladders. We have recently studied 
the possibility that equ. (74) can be justified as the leading term in an 
expansion in 1/d, where d is the number of space dimensions ^3,52 
approach appears promising, but 1/d corrections have not been worked out 
yet. Another possibility is to pursue numerical approaches. Green func- 
tion Monte Carlo calculations give ^ = 0.44, in very good agreement with 
the experimental result Several groups have performed euclidean lattice 
calculations^^'^^'^^'^*, similar to the lattice QCD calculations discussed in 
Sect. 3.4. These calculations do not suffer from a sign problem and can be 
extended to finite temperature. 

4.4. Nuclear Matter and Chiral Restoration 

Ordinary nuclei consist of roughly equal numbers of neutrons and protons. 

In order to study heavy nuclei it is useful to consider nuclear matter in pure 
QCD, i.e. ignoring the contribution from electrons as well as the Coulomb 
repulsion between protons. As discussed in Sect. 4.1 nuclear matter satu- 
rates at a density po — 0.15 fm~^. The binding energy of nuclear matter is 
B/A ~ 15 A/[eV. Numerical calculations based on realistic miclcon-nuclcon 
potentials are successful in reproducing these numbers, but we do not un- 
derstand very well why nuclear matter saturates and how the saturation 
density and the binding energy are related to more fundamental properties 
of QCD. 

We also do not know very well how to extrapolate the equation of state 
beyond the saturation density. An important question is whether we expect 
chiral symmetry to be restored as the density increases. If the density 
is small this question can be studied using the method we employed in 
Sect. 3.1. The quark condensate is given by 

The partition function of a dilute gas of protons and neutrons is 

log Z = 4 y ^ log (l + e-(^"-'')/^) . (76) 
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The quark mass dependence of the nucleon mass is related to ttN Sigma 
term S^jv = mqdniN/dmq. We get 

(77) 

The Sigma term can be extracted in pion-nucleon scattering. Using S^jv — 
45 MeV we find 

(,-,)o{l4(^)}, (78) 

which indicates that chiral condensate is significantly modified already at 
nuclear matter saturation density. 



4.5. Superfluidity 

One of the most remarkable phenomena that take place in many body 
systems is superfluidity. Superfluidity is related to an instability of the 
Fermi surface in the presence of attractive interactions between fermions. 
Let us consider fermion-fermion scattering in the simple model introduced 
in Sect. 4.2. At leading order the scattering amplitude is given by 

ral3'rSiPl,P2,P3,P'i) = Co (SafSps - SasSfj^) . (79) 

At next-to-lcading order we find the corrections shown in Fig. 11. A detailed 
discussion of the role of these corrections can be found in ^'59, 60^ r^^ie 
BCS diagram is special, because in the case of a spherical Fermi surface it 
can lead to an instability in weak coupling. The main point is that if the 
incoming momenta satisfy pi ~ —p2 then there are no kinematic restrictions 
on the loop momenta. As a consequence, all back-to-back pairs can mix 
and there is an instability even in weak coupling. 

For pi = —p2 and Ei = E2 = E the BCS diagram is given by 

^afjjS = Cq {5aiS/3S - SasSf}-^) [ JTTT^-^, -"7 

J (27r)* E + qo-eq+ idsgn{eq) 

E -qo-eq + iSsgn{eq) ' ^^^^ 

The loop integral has an infrared divergence near the Fermi surface as 
E ^ 0. The scattering amplitude is proportional to 

^a0j5 = i^a-ySpS - SasSffj) | Cq - Cq (^^) ( ) ( ' ^^^^ 
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ZS 



ZS' 



BCS 



Figure 11. Second order diagrams that contribute to particle-particle scattering. The 
three diagrams are known as ZS (zero sound), ZS' and BCS (Bardeen-Cooper-SchriefTer) 
contribution. 



where Eq is an ultraviolet cutoff. Equ. (81) can be interpreted as an effective 
energy dependent coupling that satisfies the renormalization group equation 

59,60 



j^dCo 2 {PFm\ 



with the solution 



CoiE) = 



Co{Eo 



1 



(82) 



(83) 



NCoiEo)log{Eo/E)' 

where TV = [ppm) / (2'it'^) is the density of states. Equ. (83) shows that there 
are two possible scenarios. If the initial coupling is repulsive, Cq{Eq) > 0, 
then the renormalization group evolution will drive the effective coupling 
to zero and the Fermi liquid is stable. If, on the other hand, the initial 
coupling is attractive, Co{Eo) < 0, then the effective coupling grows and 
reaches a Landau pole at 

At the Landau pole the Fermi liquid description has to break down. The 
renormalization group equation does not determine what happens at this 
point, but it seems natural to assume that the strong attractive interaction 
will lead to the formation of a fermion pair condensate. The fermion con- 
densate (c^^V'qV'/s) signals the breakdown of the U{1) symmetry and leads 
to a gap A in the single particle spectrum. 

The scale of the gap is determined by the position of the Landau pole, 
A ~ Ecrit ■ A more quantitative estimate of the gap can be obtained in the 
mean field approximation. In the path integral formulation the mean field 
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approximation is most easily introduced using the Hubbard-Stratonovich 
trick. For this purpose we first rewrite the four-fermion interaction as 

^(^t^)2 = 9l {{'4,'fa2->p^){ipa2ip) + (VV2aV^)(^a(72V)} , (85) 

where we have used the Fierz identity 25°'f^5'^P = S'^pS'^^ + {a} {aY' ^ . Note 
that the second term in equ. (85) vanishes because (cr2a) is a symmetric 
matrix. We now introduce a factor of unity into the path integral 

1 = ^/mexp(^), (86) 

where we assume that Cq < 0. We can eliminate the four-fermion term 
in the lagrangian by a shift in the integration variable A. The action is 
now quadratic in the fermion fields, but it involves a Majorana mass term 
'ipa2Aip + h.c. The Majorana mass terms can be handled using the Nambu- 
Gorkov method. We introduce the bispinor ^ = a^) and write the 

fermionic action as 

2 J {2T:f^ \ A* po + ej > 

Since the fermion action is quadratic we can integrate the fermion out and 
obtain the effective lagrangian 

L=iTV[log(Go-iG)]+^|A|2, (88) 

where G is the fermion propagator 

The diagonal and off-diagonal components of G{p) are sometimes referred to 
as normal and anomalous propagators. Note that we have not yet made any 
approximation. We have converted the fermionic path integral to a bosonic 
one, albeit with a very non-local action. The mean field approximation 
corresponds to evaluating the bosonic path integral using the saddle point 
method. Physically, this approximation means that the order parameter 
does not fluctuate. Formally, the mean field approximation can be justified 
in the large A'^ limit, where N is the number of fermion fields. The saddle 
point equation for A gives the gap equation 
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Performing the po integration we find 

2 J (2.)3 ^J^,- ^^'^ 

Since ep ~ Ep ~ ^ the integral in equ. (91) has an infrared divergence on 
the Fermi surface \p\ ^ pp. As a result, the gap equation has a non-trivial 
solution even if the coupling is arbitrarily small. The magnitude of the gap 
is A ~ Aexp(— l/(|Co|A'')) where A is a cutoff that regularizes the integral 
in equ. (91) in the ultraviolet. If we treat equ. (65) as a low energy effective 
field theory we should be able to eliminate the unphysical dependence of 
the gap on the ultraviolet cutoff, and express the gap in terms of a physical 
observable. At low density this can be achieved by observing that the gap 
equation has the same UV behavior as the Lipmann-Schwinger equation 
that determines the scattering length at zero density 



47ra 2 J {2'jTf Ep' 

Combining equs. (91) and (92) we can derive an UV finite gap equation 
that depends only on the scattering length, 

m If d?p r 1 



47ra 2 J (27r)n ^^7^:^ E, 

Solving for the A we find ^^'^^ 

A=^exp(--^). (94) 

Higher order corrections reduce the pre-exponent in this result by a factor 
(4e)^/^ ~ 2.2 Like the perturbative calculation of the energy per particle 
this result is not very useful for neutron matter, since the scattering length is 
very large. Taking higher order corrections into account, Equ. (94) suggests 
that A ~ QAQEf as pp\a\ — > oo. Surprisingly, this estimate agrees very 
well with numerical calculations The gap is also quite sensitive to the 
effective range of the interaction. Calculations based on potential models 
give gaps on the order of 2 MeV at nuclear matter density. 



4.6. Landau- Ginzburg theory 

In neutron stars there is not only pairing between neutrons but also pairing 
between protons. Since protons carry charge this implies that the material 
is not only a superfluid but also a superconductor. Superconductors have 
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many interesting properties which can be understood from the symmetries 
involved. We will consider a system of protons coupled to a U{1) gauge 
field A/j,. The order parameter $ = {e'^'^tpai^is) breaks U{1) invariance. 
Consider a gauge transformation 

Af,^Af, + di,A. (95) 

The order parameter transforms as 

$ exp(2ieA)$. (96) 

The breaking of gauge invariance is responsible for most of the unusual 
properties of superconductors 65,66 j^^^ig gggn by constructing the 

low energy effective action of a superconductor. For this purpose we write 
the order parameter in terms of its modulus and phase 

$(x) = exp(2ie0(x))$(a;). (97) 

The field corresponds to the Goldstone mode. Under a gauge transfor- 
mation (j){x) — > ^{x) + A{x). Gauge invariance restricts the form of the 
eff'ective Lagrange function 

L = ^IJ d^xF^.F^, + L,{A^ - d^(j)). (98) 

There is a large amount of information we can extract even without knowing 
the explicit form of Lg. Stability implies that A^ = dfxcj) corresponds to a 
minimum of the energy. This means that up to boundary effects the gauge 
potential is a total divergence and that the magnetic field has to vanish. 
This phenomenon is known as the Meissner effect. 

Equ. (98) also implies that a superconductor has zero resistance. The 
equations of motion relate the time dependence of the Goldstone boson 
field to the potential, 

^(a;) = -V{x). (99) 

The electric current is related to the gradient of the Goldstone boson field. 
Equ. (99) shows that the time dependence of the current is proportional 
to the gradient of the potential. In order to have a static current the 
gradient of the potential has to be constant throughout the sample, and 

the resistance is zero. 

In order to study the properties of a superconductor in more detail 
we have to specify Lg. For this purpose we assume that the system is 
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time-independent, that the spatial gradients are smaU, and that the order 
parameter is small. In this case we can write 



L, = 



(100) 

where mu and g are unknown parameters that depend on the temperature. 
Equ. (100) is known as the Landau- Ginzburg effective action. Strictly 
speaking, the assumption that the order parameter is small can only be 
justified in the vicinity of a second order phase transition. Nevertheless, 
the Landau-Ginzburg description is instructive even in the regime where 
t = (T — Tc)/Tc is not small. It is useful to decompose $ = pexp(2ie(?i>). 
For constant fields the effective potential, 

V{p) = -\m]jp^ + \gp\ (101) 

is independent of 0. The minimum is at /)§ — mj^/g and the energy den- 
sity at the minimum is given by E = —m'^/{Ag). This shows that the 
two parameters tuh and g can be related to the expectation value of $ 
and the condensation energy. We also observe that the phase transition is 
characterized by mH{Tc) = 0. 

In terms of (f) and p the Landau-Ginzburg action is given by 

L, = |d3^|_2eV(W-A)%lml,p2-i5p'-i(Vp)^|. (102) 

The equations of motion for A and p are given by 

V X i? = 4eV^ (V(?!- - Aj , (103) 
VV = -m%p'^ + gp^ + 4eV (V0 - Aj . (104) 

Equ. (103) implies that V^-B = —Ae^p^B. This means that an exter- 
nal magnetic field B decays over a characteristic distance A = \/{2ep). 

Equ. (104) gives V^p = —rnjjp + As a consequence, variations in the 

order parameter relax over a length scale given by ^ = l/mjj. The two 
parameters A and ^ are known as the penetration depth and the coherence 
length. 

The relative size of A and ^ has important consequences for the prop- 
erties of superconductors. In a type II superconductor ^ < A. In this case 
magnetic flux can penetrate the system in the form of vortex lines. At the 
core of a vortex the order parameter vanishes, p = 0. In a type II material 
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the core is much smaller than the region over which the magnetic field goes 
to zero. The magnetic flux is given by 



and quantized in units of Tr/i/e. In a type II superconductor magnetic 
vortices repel each other and form a regular lattice known as the Abrikosov 
lattice. In a type I material, on the other hand, vortices are not stable 
and magnetic fields can only penetrate the sample if superconductivity is 
destroyed. 

5. QCD at high density 
5.1. Color superconductivity 

In Sect. 4.1 we introduced a few simple arguments concerning the phase 

diagram of QCD in the /i — T plane. These arguments are summarized in 
Fig. 7. The basic idea is that large baryon density is just like high temper- 
ature: there is a large scale in the problem, the eflFective coupling is weak, 
and the system is described, to a good approximation, as a weakly inter- 
acting quark liquid. We expect, in particular, that quarks are deconfined 
and that chiral symmetry is restored. 

We also showed, however, that systems at finite density, as exemplified 
by nuclear matter, have a very rich phase diagram. We saw, in particular, 
that the BCS instability will lead to pair condensation whenever there is an 
attractive fermion-fermion interaction, even if the interaction is weak. At 
very large density, the attraction is provided by one-gluon exchange between 
quarks in a color anti-symmetric 3 state. High density quark matter is 
therefore expected to behave as a color superconductor 67,68,69, 70_ 

Color superconductivity is described by a pair condensate of the form 



Here, C is the charge conjugation matrix, and Td,^c,tf are Dirac, color, 
and flavor matrices. Except in the case of only two colors, the order pa- 
rameter cannot be a color singlet. Color superconductivity is therefore 
characterized by the breakdown of color gauge invariance. This statement 
has to be interpreted in the sense of Sect. 4.6. Gluons acquire a mass due 
to the (Meissner- Anderson) Higgs mechanism. 

A rough estimate of the critical density for the transition from chiral 
symmetry breaking to color superconductivity, the superconducting gap 




(105) 



(106) 
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Figure 12. First revision of the phase diagram of hadronic matter. This figure shows 
the phase diagram of strongly interacting matter obtained from a mean field treatment 
of chiral symmetry breaking and color superconductivity in QCD with two flavors. 



and the transition temperature is provided by schematic four-fcrmion mod- 
els '^^^'''2. Typical models are based on the instanton interaction 

C = Gi{{^T-^f + {4,^^T-i^f], (107) 

or a schematic one-gluon exchange interaction 



C = GoGE[i^l^.^i^j . (108) 

Here t~ = (r, i) is an isospin matrix and A" are the color Gell-Mann ma- 
trices. The strength of the four-fermion interaction is typically tuned to 
reproduce the magnitude of the chiral condensate and the pion decay con- 
stant at zero temperature and density. In the mean field approximation 
the effective quark mass associated with chiral symmetry breaking is deter- 
mined by a gap equation of the type 

= Gm t ,f^, , (1 - np{E,)) , (109) 

where Gm is the effective coupling in the quark-anti-quark channel, A is 
a cutoff, and np{E) is the Fermi distribution. Both the instanton inter- 
action and the one-gluon exchange interaction are attractive in the color 
anti-triplet scalar diquark channel e'^^'^{tp''G^5tp'^). A pure one-gluon ex- 
change interaction leads to a degeneracy between scalar and pseudoscalar 
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Figure 13. Location of the (pseudo) critical line and tri-critical point (box) as mea- 
sured in different simulations, de Forcrand and Philipsen (2002) [FP], Fodor and Katz 
(2002) [FK], AUton et al (2002) [All], D'EUia and Lombardo (2002) [EL]. Figure from 
de Forcrand and Philipsen (2003). 



diquark condensation, but instantons are repulsive in the pseudoscalar di- 
quark channel. The gap equation in the scalar diquark channel is 

A.^r^ ^ (110) 

where we have neglected terms that do not have a singularity on the Fermi 
surface |p| = pf- In the case of a four-fermion interaction with the quantum 
numbers of one-gluon exchange Go = Gm/{Nc— 1). The same result holds 
for instanton effects. In order to determine the correct ground state we 
have to compare the condensation energy in the chiral symmetry broken 
and diquark condensed phases. We have E ~ /^Mq in the {qq) condensed 
phase and E ~ p|,A^/(27r^) in the {qq) condensed phase. 

At zero temperature and density both equs. (109) and (110) only have 
non-trivial solutions if the coupling exceeds a critical value. Since Gm > 
Gd we have Mq > A and the energetically preferred solution corresponds 
to chiral symmetry breaking. If the density increases Pauli-Blocking in 
equ. (109) becomes important and the effective quark mass decreases. The 
diquark gap equation behaves very differently. Equ. (110) has an infrared 
singularity on the Fermi surface, p = pp, and this singularity is multiplied 
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by a finite density of states, A'' = p|./(27r)^. As a consequence, there is 
a non-trivial solution even if the coupling is weak. The gap grows with 
density until the Fermi momentum becomes on the order of the cutoff. 
For realistic values of the parameters we find a first order transition for 
remarkably small values of the quark chemical potential, /ig ~ 300 McV. 
The gap in the diquark condensed phase is A ~ 100 MeV and the critical 
temperature is To ~ 50 MeV. 

In the same model the finite temperature phase transition at zero baryon 
density is found to be of second order. This result is in agreement with uni- 
versality arguments and lattice results. If the transition at finite density 
and zero temperature is indeed of first order then the first order transition 
at zero baryon density has to end in a tri-critical point 'i'4>75,76 rpj^^ 
critical point is quite remarkable, because it remains a true critical point 
even if the quark masses are not zero. A non-zero quark mass turns the 
second order T ^ transition into a smooth crossover, but the first order 
/i ^ transition persists. It is hard to predict where exactly the tri-critical 
point is located in the phase diagram. Recent lattice calculations suggest 
that the tri-critical point is sufiiciently close to the finite temperature axis 
so that its location can be determined on the lattice, sec Fig. 13. It may 
also be possible to locate the critical point experimentally. Heavy ion colli- 
sions at relativistic energies produce matter under the right conditions and 
experimental signatures of the tri-critical point have been suggested in 

A schematic phase diagram is shown in Fig. 12. We should emphasize 
that this phase diagram is based on simplified models and that there is 
no proof that the transition from nuclear matter to quark matter along 
the T = line occurs via a single first order transition. Chiral symmetry 
breaking and color superconductivity represent two competing forms of or- 
der, and it seems unlikely that the two phases are separated by a second 
order transition. However, since color superconductivity modifies the spec- 
trum near the Fermi surface, whereas chiral symmetry breaking operates 
near the surface of the Dirac sea, it is not clear that the two phases can- 
not coexist. Indeed, there are models in which a phase coexistence region 
appears 

5.2. Phase structure in weak coupling 
5.2.1. QCD with two flavors 

In this section we shall discuss how to use weak coupling methods in order 
to explore the phases of dense quark matter. We begin with what is usu- 
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Figure 14. Fig. a) shows the leading contribution to the Dyson-Schwinger (gap) equa- 
tion in QCD at finite density. The open square denotes an anomalous self energy (gap) 
insertion and the solid square is a gluon self energy insertion. Figs, b) and c) show quark 
self energy insertions and vertex corrections. 

ally considered to be the simplest case, quark matter with two degenerate 
flavors, up and down. Renormalization group arguments suggest ''^.so^ g^^^ 
explicit calculations confirm 8i,82^ ^]^g^^ whenever possible quark pairs con- 
dense in an s-wave state. This means that the spin wave function of the pair 
is anti-symmetric. Since the color wave function is also anti-symmetric, the 
Pauli principle requires the flavor wave function to be anti-symmetric too. 
This essentially determines the structure of the order parameter '^^^'^^ 

^ (gafcc^bCTsraV-^). (Ill) 

This order parameter breaks the color SU (3) — > SU (2) and leads to a gap 
for up and down quarks with two out of the three colors. Chiral and isospin 
symmetry remain unbroken. 

We can calculate the magnitude of the gap and the condensation energy 
using weak coupling methods. In weak coupling the gap is determined by 
ladder diagrams with the one gluon exchange interaction. These diagrams 
can be summed using the gap equation 83,84,85,86,87 

A(p4) ^ [ dq^ [ dcos0 ( (112) 



1 I 1 



A(g4) 



I ' cose + F/ (2(1^) J ^ql + Aiq^r' 

Here, A(p4) is the frequency dependent gap, g is the QCD coupling constant 
and G and F are the self energies of magnetic and electric gluons. This gap 
equation is very similar to the BCS gap equation equ. (110) obtained in four- 
fermion models. The terms in the curly brackets arise from the magnetic 
and electric components of the gluon propagator. The numerators are the 
on-shell matrix elements Mu^qq = [uhiPihiflUh{p3)][uh{p2hi,QUhiP4:)] for 
the scattering of back-to-back fermions on the Fermi surface. The scattering 
angle is cos 9 = pi -p^. In the case of a spin zero order parameter, the helicity 
h of all fermions is the same, see for more detail. 
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The main difference between equ. (112) and the BCS gap equation (110) 
is that because the ghion is niasslcss, the gap equation contains a colhnear 
divergence for cos^ ~ 1. In a dense medium the coUinear divergence is 
regularized by the gluon self energy. For g — > and to leading order in 
perturbation theory we have 

F = 2m^ G=-m''%, (113) 

with = Njg^ jj? / {Att"^). In the electric part, m|, = 2m^ is the familiar 
Debyc screening mass. In the magnetic part, there is no screening of static 
modes, but non-static modes arc modes arc dynamically screened due to 
Landau damping. This is completely analogous to the situation at finite 
temperature see Sect. 3.5. 

For small energies dynamic screening of magnetic modes is much weaker 
than Debye screening of electric modes. As a consequence, perturbative 
color superconductivity is dominated by magnetic gluon exchanges. Using 
equ. (113) we can perform the angular integral in equ. (112) and find 

with b = 256n^{2/Nf)^/'^g~^. We can now see why it was important to keep 
the frequency dependence of the gap. Because the coUinear divergence is 
regulated by dynamic screening, the gap equation depends on even if the 
frequency is small. We can also see that the gap scales as exp{—c/g). The 
coUinear divergence leads to a gap equation with a double-log behavior. 
Qualitatively 

^-il^^h(f)]'' ^'''^ 

from which we conclude that A ~ exp{—c/g). Equ. (115) is not suflaciently 
accurate to determine the correct value of the constant c. A more detailed 
analysis shows that the gap on the Fermi surface is given by 

Ao ~ 512n^{2/Nff/%'ng-^ exp • (US) 

The factor b' is related to non- Fermi liquid effects, see 89,90,9i_ jj^ pgp^;ur- 
bation theory b' = exp{-{n'^ + 4){Nc - 1)/16) 87,92,93 r^^ie condensation 
energy is given by 
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where Nd = 4 is the number of condensed species. In the mean field 
approximation the critical temperature is Tc/Aq = e* j-K ~ 0.56, as in 
standard BCS theory Fluctuations of the gauge field drive the transition 
first order Surprisingly, gauge field fluctuations increase the critical 

temperature as compared to the BCS result 

For chemical potentials /x < 1 GeV, the coupling constant is not small 
and the applicability of perturbation theory is in doubt. If we ignore this 
problem and extrapolate the perturbative calculation to densities p ~ 5/9o 
we find gaps on the order of lO's of MeV. However, perturbative estimates 
also show that instanton effects cannot be neglected for < 1 GeV, and 
that instantons increase the gap 

We note that the 2SC phase defined by equ. (Ill) has two gapless 
fermions and an unbroken SU{2) gauge group. The gapless fermions are 
singlets under the unbroken SU (2). As a consequence, we expect the SU (2) 
gauge group to become non-pcrturbativc. An estimate of the SU (2) con- 
finement scale was given in We also note that even though the Copper 
pairs carry electric charge the t/(l) of electromagnetism is not broken. The 
generator of this symmetry is a linear combination of the original electric 
charge operator and the diagonal color charges. Under this symmetry the 
gapless fermions carry the charges of the proton and neutron. Possible 
pairing between the gapless fermions was discussed in '^i>9*. 

5.2.2. QCD with three flavors: Color-Flavor-Locking 

If quark matter is formed at densities several times nuc;lear matter density 
we expect the quark chemical potential to be larger than the strange quark 
mass. We therefore have to determine the structure of the superfluid order 
parameter for three quark flavors. We begin with the idealized situation of 
three degenerate flavors. From the arguments given in the last section we 
expect the order parameter to be color-flavor matrix of the form 

<1>^/ = (CC^75^,'). (118) 

The structure of this matrix can be determined by extremizing the grand 
canonical potential. We flnd ^^-^oo 

A«/ = A^(<5«5) - 6^6]) + As{S16'^ + Sl6]), (119) 

which describes the color-flavor locked (CFL) phase proposed by Alford, 
Rajagopal, and Wilczck In the weak coupling limit Ag <C and 
Aa = 2~i/^Ao where Aq is the gap in the 2SC phase, equ. (116) In 
the CFL phase both color and flavor symmetry are completely broken. 
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Figure 15. Conjectured phase diagram of Nf = 3 hadronic matter in the limit of exact 
flavor symmetry. 

There are eight combinations of color and flavor symmetries that generate 
unbroken global symmetries. The unbroken symmetries are 

- {UT'U.^^lr rn,. - iU*rU.,^Hj, (120) 

for U G SU{3)v- The symmetry breaking pattern is 

SU{3)l X SU{3)h X U{l)v SU{3)v (121) 

We observe that color-flavor- locking implies that chiral symmetry is broken. 
The mechanism for chiral symmetry breaking is quite unusual. The primary 
order parameter (tpl^CAf-'tlj'lj) = —{tp'j^^CA'^j^p^j) involves no couphng 
between left and right handed fermions. In the CFL phase both left and 
right handed flavor are locked to color, and because of the vectorial coupling 
of the gluon left handed flavor is effectively locked to right handed flavor. 
Chiral symmetry breaking also implies that {^pip) has a non-zero expectation 
value In the CFL phase {ip'ip)'^ ^ {{ipip)'^). Another measure of chiral 
symmetry breaking is provided by the pion decay constant. We will see 
that in the weak coupling limit is proportional to the density of states 
on the Fermi surface. 

The symmetry breaking pattern SU{3)l x SU{3)ii SU{3)v in the 
CFL phase is identical to the symmetry breaking pattern in QCD at low 
density. The spectrum of excitations in the color-flavor-locked (CFL) phase 
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also looks remarkably like the spectrum of QCD at low density The 
excitations can be classified according to their quantum numbers under 
the unbroken SU{3), and by their electric charge. The modified charge 
operator that generates a true symmetry of the CFL phase is given by 
a linear combination of the original charge operator Qem, and the color 
hypercharge operator Q = diag(— 2/3, —2/3, 1/3). Also, baryon number is 
only broken modulo 2/3, which means that one can still distinguish baryons 
from mesons. We find that the CFL phase contains an octet of Goldstone 
bosons associated with chiral symmetry breaking, an octet of vector mesons, 
an octet and a singlet of baryons, and a singlet Goldstone boson related to 
superfluidity. All of these states have integer charges. 

With the exception of the U{1) Goldstone boson, these states exactly 
match the quantum numbers of the lowest lying multiplets in QCD at low 
density. In addition to that, the presence of the U{1) Goldstone boson 
can also be understood. The U{1) order parameter is {{uds){uds)) . This 
order parameter has the quantum numbers of a 0"*" (AA) pair condensate. 
In Nf = 3 QCD, this is the most symmetric two nucleon channel, and a 
very likely candidate for superfluidity in nuclear matter at low to moderate 
density. We conclude that in QCD with three degenerate light flavors, 
there is no fundamental difference between the high and low density phases. 
This implies that a low density hyper-nuclear phase and the high density 
quark phase might be continuously connected, without an intervening phase 
transition. A conjectured phase diagram is shown in Fig. 15. 

5.3. The role of the strange quark mass 

At baryon densities relevant to astrophysical objects distortions of the pure 
CFL state due to non-zero quark masses cannot be neglected i04,i05_ -pj^g 
most important efii'ect of a non-zero strange quark mass is that the light 
and strange quark Fermi momenta will no longer be equal. When the 
mismatch is much smaller than the gap one calculates assuming degenerate 
quarks, we might expect that it has very little consequence, since at this 
level the original particle and hole states near the Fermi surface are mixed 
up anyway. On the other hand, when the mismatch is much larger than 
the nominal gap,we might expect that the ordering one would obtain for 
degenerate quarks is disrupted, and that to a first approximation one can 
treat the light and heavy quark dynamics separately. 

This argument is qualitatively right, but the correct picture turns out 
to be much more complicated, and much more interesting. If the strange 
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Figure 16. Hierarchy of effective field theories in the CFL phase. 



quark mass is taken into account microscopic calculations based on the 
Dyson-Schwinger equation become much more complicated, because there 
are many more gap parameters, and maintaining electric neutrality and 
color gauge invariance is difficult I06,i07,i08^ However, since chiral symmetry 
is broken in the CFL phase we know that the dependence on the quark 
masses is constrained by chiral symmetry. It is therefore natural to study 
the problem using effective field theories. In practice we will employ a two- 
step procedure, see Fig. 16. In the first step we match the microscopic 
theory, QCD, to an effective field theory of quasi-particles and holes in the 
vicinity of the Fermi surface. In the second step we match this theory to 
an effective chiral theory for the CFL phase. 

5.3.1. High density effective theory 

The QCD Lagrangian in the presence of a chemical potential is given by 



where Z?^ = 9^ + igA^ is the covariant derivative, M is the mass ma- 
trix and fi is the baryon chemical potential. If the baryon density is very 
large perturbative QCD calculations can be further simplified. The main 
observation is that the relevant degrees of freedom are particle and hole 
excitations in the vicinity of the Fermi surface. We shall describe these 
excitations in terms of the field ip+ {v, x) , where v is the Fermi velocity. At 
tree level, the quark field "0 can be decomposed as ip = ^+ -I- "0- where 



£ = ill/) + M7o) ^ - i'LM^R - ^rMHl - -^G%Gl, 



(122) 
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Figure 17. Mass terms in the high density effective theory. The first diagram shows 
a O(MM^) term that arises from integrating out the tfi— field in the QCD lagrangian. 
The second diagram shows a O(M^) four-fermion operator which arises from integrating 
out ip— and hard gluon exchanges. 

ijj± = i(l ± (3 • v)iIj. Note that (1 ± c? • ?))/2 is a projector on states with 
positive/negative energy. To leading order in 1/pf we can eliminate the 
field ip- using its equation of motion. The lagrangian for the 1/'+ field is 
given by W9awmi 

C = ^jliiv ■ D)^+ - ^ (t^^'Ci/;^^' {5a^5bJ - 5aj5u) + h.c.) + . . . , (123) 

with — (1, v) and . . . and a, 6, . . . denote flavor and color indices. The 
magnitude of the gap A is determined order by order in perturbation theory 
from the requirement that the thermodynamic potential is stationary with 
respect to A. With the gap term included the perturbative expansion is 
well defined. 

The effective theory contains an infinite set of operators that have ad- 
ditional derivatives or more powers of "0+ . These operators are suppressed 
by inverse powers of the Fermi momentum. Here, we will only consider op- 
erators that contain quark masses. To leading order in 1/pf there is only 
one operator in the high density effective theory 

C = (vl+MMtV'i+ -f V;^^+MtMV'i^+) . (124) 

This term arises from expanding the kinetic energy of a massive fermion 
around p ~ pp. At 0(l/p|,) we find four-fermion operators that contain 
two powers of the quark mass. The coefficients of these operators are ob- 
tained by computing chirality violating quark-quark scattering amplitudes 
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for quasi-particles near the Fermi surface see Fig. 17. At leading order 
in 1/pf these amphtudcs are independent of the scattering angle and can 
be represented as local four-fcrmion operators 

(125) 

There are two additional terms with (L R) and (M ^ M^). We have 
introduced the CFL eigenstates tp^ defined by ip"; — 'ip^{\^)ai/^/^, A = 
0, . . . , 8. The tensor F is defined by 

J.ABCD ^ 1 [\^M{\^f\BM{\^f] 

- \tv [A^M(A^)^] Tr [A^M(A^)^] |. (126) 
The explicit expression for F is given in but we will not need here. 

5.3.2. CFL chiral theory 

For excitation energies smaller than the gap the only relevant degrees of 
freedom are the Goldstone modes associated with the breaking of chiral 
symmetry and baryon number, see Fig. 16. Since the pattern of chiral 
symmetry breaking is identical to the one at T = = the effective 
lagrangian has the same structure as chiral perturbation theory. The main 
difference is that Lorentz-invariance is broken and only rotational invariance 
is a good symmetry. The effective lagrangian for the Goldstone modes is 

113 

£e// = y Tr [VoSVoS^ - vldiY.diY.'^] + [BTr(MI]t) + h.c] 
+ [AiTr(MS^)Tr(MEt) + y42Tr(MStMSt) 

+yl3Tr(MSl')Tr(Ml'S) + h.c] + .... (127) 

Here S = exp{i(l)°- X°- / Jt^) is the chiral field, /tt is the pion decay constant 
and M is a complex mass matrix. The chiral field and the mass matrix 
transform as E ^ LTiR^ and M LMR^ under chiral transformations 
{L, R) e SU (3)l X SU {3)r. We have suppressed the singlet fields associated 
with the breaking of the exact U{l)v and approximate U{1)a symmetries. 

At low density the coefficients /,r, B,Ai,... are non-perturbative quan- 
tities that have to extracted from experiment or measured on the lattice. 
At large density, on the other hand, the chiral coefficients can be calculated 
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in perturbative QCD. The leading order terms are 



114 



, 21-81og(2) 1 

U- [2^)' ""--s- ^^^^^ 

Mass terms arc determined by the operators studied in the previous section. 
We observe that both equ. (124) and (125) are quadratic in M. This implies 
that B = in perturbative QCD. B receives non-perturbative contributions 
from instantons, but these effects are small if the density is large 

We observe that Xl = MMt/(2pF) and Xr = M^M/{2pf) in 
equ. (124) act as effective chemical potentials for left and right-handed 
fermions, respectively. Formally, the effective lagrangian has an SU{3)l x 
SU{3) u gauge symmetry under which X^ ji transform as the temporal com- 
ponents of non-abelian gauge fields. We can implement this approximate 
gauge symmetry in the CFL chiral theory by promoting time derivatives to 
covariant derivatives 

/MM''\ /M''M\ 

v.E.a,E + i(— )e-.e(_). (m, 

The four-fermion operator in equ. (125) contributes to the coefficients 
Ai.We find "4,112 

Ai = -A2 = ^3 = 0. (130) 

We can now summarize the structure of the chiral expansion in the CFL 
phase. The effective lagrangian has the form 

Loop graphs in the effective theory are suppressed by powers of d/ {4:TTf^). 
Since the pion decay constant scales as /,r ^ Goldstone boson loops 

are suppressed compared to higher order contact terms. We also note that 
the quark mass expansion is controlled by m^/ (ppA), as expected from the 
arguments presented in Sect. 5.3. 



5.3.3. Kaon condensation 

Using the chiral effective lagrangian we can now determine the dependence 
of the order parameter on the quark masses. We will focus on the physically 
relevant case rrig > to„ = m^- Because the main expansion parameter is 
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Figure 18. This figure shows the phase structure of CFL matter as a function of the 
strange quark mass rrts and the lepton chemical potential ^q, from Kaplan and Reddy 
(2001). 



mg/(pi?A) increasing the quark mass is roughly equivalent to lowering the 
density. The effective potential for the order parameter is 



£2 

Veff = yTr [Xl^XrT.^ + Ai (Tr(MI]^))^ - Tr(MEtMS^) 



If the strange quark mass is small then the minimum of the effective po- 
tential is E = 1. However, when the strange quark mass exceeds a certain 
critical value it becomes favorable to rotate the order parameter in the 
kaon direction. The physical reason is that the system tries to reduce its 
strangeness content by forming a kaon condensate. Consider the ansatz 
E = exp(zQ;A4). The vacuum energy is 

V{a) = [ ^ ( ) sin(a)2 + (mS,)2(cos(a) - 1) ) , (133) 




where (m^f)^ ~ {4:Ai / f^)m{m + nis). Minimizing the vacuum energy we 
obtain a = if /is < and cos(q!) = {i^^kY I A if Ms > Here, 
we have defined fig — ra1/{2pp). Using the perturbative result for Ai the 
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20 40 60 80 100 

mll{2pF) [MeV] 

Figure 19. This figure shows the fermion spectrum in the CFL phase. For ms = 
there are eight fermions with gap A and one fermion with gap 2 A (not shown). Without 
kaon condensation gapless fermion modes appear at /i^ = A (dashed hnes). With kaon 
condensation gapless modes appear at fis = 4A/3. 

critical strange quark mass is 

m,|_, = 3.03-my^A2/3. (134) 

Using A ~ 50 MeV we get rrigicrit) ~ 70 MeV. This result suggests that 
strange quark matter at densities p ^ {h — 10)/9o is in a kaon condensed 
phase. The kaon condensate breaks SU{2)j x U{1)y to U{1)q. The phase 
structure as a function of the strange quark mass and non-zero lepton chem- 
ical potentials was studied by Kaplan and Reddy see Fig. 18. We ob- 
serve that if the lepton chemical potential is non-zero charged kaon and 
pion condensates are also possible. It was also shown that there is a range 
of light quark masses in which simultaneous kaon and eta condensation 
takes place 

6. Fermions in the CFL phase 

So far we have only studied Goldstone modes in the CFL phase. However, as 
the strange quark mass is increased it is possible that some of the fermion 
modes become light or even gapless In order to study this question 
we have to include fermions in the effective field theory. The effective 



Februciry 2, 2008 3:55 Proceedings Trim Size: 9in x 6in 



hugs'proc 



50 

lagrangian for fermions in the CFL phase is i^^'^^o 

£ = Tr (nUv^Di^N) - DTr {N^v^^'y^ {A^, N}) - FTr {Nh^^js N]) 

+ ^ { (Xr (NlNl) - [Tr (iV^)]') - {L ^ R) + h.c.} . (135) 

Nl,r are left and right handed baryon fields in the adjoint representation of 
flavor SU{3). The baryon fields originate from quark-hadron complemen- 
tarity as explained in Sect. 5.2.2. We can think of N as being composed 
of a quark and a diquark field, Nl ~ qL{qLqL)- The covariant derivative 
of the nucleon field is given by f^AT = d^N + i[V^, N]. The vector and 
axial- vector currents are 

where ^ is defined by = E. It follows that ^ transforms as ^ ^ L^U {xY = 
U{x)^R^ with U{x) e SU{3)v. For pure SU{3) flavor transformations 

L = R = V wc have U{x) = V. F and D are low energy constants 
that determine the baryon axial coupling. In perturbative QCD we find 
D = F = 1/2. 

Mass terms can be introduced as in Sect. 5.3.2. The {Xl, Xr) covariant 
derivative of the nucleon field is 

DoN = doN + i[To,N\, (137) 
To = {do + iXn) + (ao + iXl) , 

where Xl = MMt/(2pF) and Xr = M'^M/{2pf) as before. The spectrum 
of fcrmion is shown in Fig. 19. A gaplcss fcrmion mode appears at /i^ = 
4A/3. In the vicinity of this point the homogeneous CFL phase becomes 
unstable 121,122 eff'ective field theory this manifests itself as an 

instability with respect to the generation of a non-zero current 123,124_ Ytom 
the effective lagrangian equ. (127) wc see that a meson current has energy 
£ ~ f^j^. This is not the end of the story, however, because a meson 
current also modifies the fermion dispersion relation. The energy of the 
lowest mode in the background of a hypercharge current Jk is given by 

;2 3 1 ^ 
'^'"'^+2A"4^^"4'^'-^''^' ^^^^^ 
where / is the momentum relative to the Fermi surface. We observe that 
the current lowers the energy of the fermions on part of the Fermi sur- 
face. When these states become gapless the total energy is lowered and 
the system can become unstable. The new ground state is a p-wave meson 
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highT 




high energy large density 

color glas color superconductivity 

Figure 20. The many facets of QCD. 

condensate in which the non-zero meson current is balanced by a back- 
flow of gapless fermions. At even larger values of this state may evolve 
into an inhomogeneous superconductor of the type considered by Larkin, 
Ovchinnikov, Fulde and Ferrell 125,126^ 127^ 



7. Outlook 

Figure 20 indicates that there are many interesting connections that we have 
not been able to explore in these lectures. There has been a lot of progress 
in connecting the large Nc world to a theory of strings, and this connection 
also sheds some light on the behavior of a strongly coupled QCD plasma. 
The transport properties of the strongly coupled plasma are probably quite 
similar to the transport behavior of the strongly correlated neutron fluid, 
and this system is related to cold trapped fermionic atoms near a Feshbach 
resonance. A lot of progress has been made in understanding hot and dense 
pre-equilibrium states, and these states share some of the properties of 
equilibrium phases. Many more surprising connections are likely to energy 
in the future. 

Acknowledgments: I would like to thank the organizers of the HUGS 
summer school, especially Jose Goity, for their hospitality. The original lec- 
tures contained a summary of the experimental work at RHIC and CERN. 
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This material is not included in the write-up, but the slides are still available 
on my website. The second half of the lectures is an abridged and updated 
version of the NPSS lectures Fig. 20 was inspired by R. Brower. This 
work was supported in part by a US DOE grant DE-FG-88ER40388. 
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